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ABSTRACT 

We report new Hersche^ observations of 25 z — 4.8 extremely luminous optically selected active galactic 
nuclei (AGNs). Five of the sources have extremely large star forming (SF) luminosities, Lsf, corresponding 
to SF rates (SFRs) of 2800-5600 M ra yr~' assuming a Salpeter IMF. The remaining sources have only upper 
limits on their SFRs but stacking their Herschel images results in a mean SFR of 700 ± 150 Mq yr _1 . The 
higher SFRs in our sample are comparable to the highest observed values so far, at any redshift. Our sample 
does not contain obscured AGNs, which enables us to investigate several evolutionary scenarios connecting 
super-massive black holes and SF activity in the early universe. The most probable scenario is that we are 
witnessing the peak of SF activity in some sources and the beginning of the post-starburst decline in others. We 
suggest that all 25 sources, which are at their peak AGN activity, are in large mergers. AGN feedback may be 
responsible for diminishing the SF activity in 20 of them but is not operating efficiently in 5 others. 
Subject headings: galaxies: active — galaxies: star formation — quasars: general 



1. INTRODUCTION 

The growth and evolution of super-massive black holes 
(SMBHs) are closely connected to the formation and evolu- 
tion of their host galaxies. While SMBHs grow through accre- 
tion of matter from their surroundings during an active galac- 
tic nucleus (AGN) phase, their host galaxies grow by star for- 
mation (SF). The bolometric luminosity of the AGN (Lagn) 
is found to be related to the SF luminosity (Lsf; defined as 
the integrated luminosity due to SF between 8-1000 /xm) 
in its host galaxy. In AGN-dominated sources (those with 
Lagn > ^SfX there is evidence for a simple power-law rela- 
tionship between the two luminosities that can be expressed 
as L SF - 10 43 ('£Ar.M/('10 43 er gs- 1 'l'l - 7 (iNetzerl 120091). This 
relationship is fou nd in tvpe-I dNetzer et al.l l2007b; Lu tz et alj 
2008) and type-II dNetzei 
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sources. 



Recent studies prov ide mo re insight into the gen eral Lagn- 
Lsf dependence. Sh ao et all (1201 Oh andlHatziminaoglo u et alJ 
(2010) measured Lsf for the hosts of known AGNs in 
the GOODS-N and HerMES fields, respectively. Both 
studies find indications that in SF-dominated sources (i.e. 
Lsf ^Lagn) there is no c lear correla t ion b etween the two 
luminosities. Furthermore, IShao et all (120101) found that for 
SF-dominated sources, in a given redshift bin, Lsf is roughly 
constant with increasing Lagn- This behavior can be inter- 
preted as an indication for different stages of evolution, where 
in sources with Lsf > Lagn the SMBH has not yet reached 
its highest accretion rate phase. These Herschel-based stud- 
ies included mostly low Lagn sources and very few sources at 
z >2.5. At higher redshifts and higher Lagn there are very few 
known SF-dominated sources although a handful of sources 
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with Lsf— L actN have been found mainly by sub-mm obser- 
vations (e.g., Isaak et al. 2002| iPriddev et al]|2003t IWu et alj 
2009; Leipski et al. 2010). 

G alaxy evolution scenario s suggest two modes of SF (see 
e.g.,|Rodighiero et al1 l2011l) . The steadier process of secular 
evolution is common in isolated disk galaxies and can reach 
SF rates (SFRs) of ~ 400 M e yr _1 at high redshift. A less 
common process with SFR that can exceed ^1000 M Q yr _1 , 
is associated with m ergers between two l arge gas-rich galax- 
ies (iDi Matteo et al] 12005b iGuyon etal] 120061: Isljacki et alJ 
201 1; Valiante et al. 201 1). Both processes result in cold gas 



inflow into the center of the system which can trigger AGN 
activity. Numerical simulations and semi-analytic models of 
mergers suggest that in such events, the fastest SMBH growth 
phase succeeds the peak of SF activity by several hundred 
Myr and is likely to take place when the SMBH is obscured 
(e.g. JHopkinretani2006l IDi Matteo et alJl2008l) . If correct, 
it means that the most luminous SF phase of the merger pre- 
cedes the most luminous AGN phase. Observationally this is 
supported by the fact that sub-mm galaxies (SMGs) have high 
Lsf, however, they often exhibit little or no AGN activity. 

SF and AGN activity may also be related through AGN 
feedback. This process can diminish or even terminate SF 
and SMBH accretion through fast winds and intense AGN 
radiation (e.g., Pi Matteo et al. 2005; Spr ingel et al.1 120051; 
Siiack fet alj I2007I) . The power of AGN feedback depends 
on Lagn and can be important in the final stages of both SF 
modes. 

In this Letter we report new Herschel observations of 25 
optically selected AGNs f rom our z — 4.8, flux-limited sample 
(Trakhtenbrot et al. 201 1, hereafter Til). Our results provide 
evidence for extreme SFRs in 20% of the sources, indicating 
merger-driven SF activity. In section|2]we describe the sample 
and the observations and explain our method for deriving Lsf 
and SFR. Section[3]compares our findings with earlier results 
and discusses the way they improve the understanding of the 
relations between AGN activity and SF in the early universe. 
Throughout this Letter we assume Hq = 70 kms~' Mpc -1 , 
D. M = 0.3 and D. A = 0.7. 



TABLE 1 

L S F AND SFRS FOR THE FIVE INDIVIDUALLY DETECTED AND STACKED SOURCES 



Object name 


Redshift 


Lagn" 


/250 


/350 


/500 


Z.SF* 


SFR'' 






(10 47 ergs-': 


1 mjy 


mJy 


mJy 


(10 l3 L a ) 


(M yi- 1 ) 


J0331-0741 


4.73 


1.23 


27.6 ± 6.2 


29.1 ±6.7 


19.1 ±7.0 


7 1Q+D-38 


377 1±$! 


J0807+1328 


4.88 


1.17 


13.5 ±5.9 


22.7 ± 6.5 


21.3 ±7.0 


1 65 +0 ' 45 


2860±™ 


J1341+0141 


4.69 


1.82 


34.5 ± 6.4 


44.0 ±7.1 


37.0 ±7.4 


3 26 +a39 
-'■ zu -0.28 


5645+*^ 


J1619+1238 


4.81 


0.95 


39.4 ± 6.5 


32.9 ±6.8 


21.2 ±7.0 


3 00 +0 - 20 


5184j*« 


J2225-0014 


4.89 


1.70 


23.3 ±6.1 


26.2 ± 6.6 


18.6 ±6.9 


2.1S±l 


377lii 


Stacked source 


4.75 


0.78 


4 53+ 1 - 4 
H - J -1.02 


6 22 +Le,/L 


3-65±i;« 


0.32j£_* 


697-- 



" Values taken from Tl 1 except J1619+1238 
* Assuming that there is no confusing source 
c Assuming a Salpeter IMF 



for which Lagn was calculated as described in Tl 1. 
within 10" from the AGN 



2. OBSERVATIONS AND DATA ANALYSIS 
2. 1 . Basic Measurements 

The 25 sources presented here are part of our sample of 40 
Z ~ 4.8 luminous, unobscured, AGNs selected from the Sloan 
Digital Sky Survey (SDSS) and described in detail in Til. 
The sources were selected from the SDSS/DR6 database, re- 
quiring 4.65 < z < 4.92 and f x (1450 A)> 6 x lO^^ergs" 1 
A -1 , to ensure sufficient H-band brightness. Since the sources 
were optically selected they are not biased with regards to SF 
properties. The redshift of 4.8 was chosen to allow the mea- 
surement of the Mg II A 2798 emission line and the AGN con- 
tinuum luminosity at rest-frame 3000 A, using H-band spec- 
troscopy. These allow a reliable measurement of Mbh and 
normalized accretion rate, L/Ledd (McLure & Dunlop 2004). 
Based on this sample Til showed that the highest lumi- 
nosity AGNs at z — 4.8 have, on average, lower Mbh and 
higher L/Lg^d than the highest luminosity AG Ns in the corre- 
sponding s amples of sources at z ~2-3.5 from lShemmer et all 
(2004) and. Netzeretal.l(l2007al) . Most of the z ~ 4.8 SMBHs 
seem to be at the end of their first continuous growth phase 
that started at z ^ 10 and are on their way to become the most 
massive (> 10 10 M @ ) black holes in the Universe. An addi- 
tional object, J1619+1238, was selected as part of the z — 4.8 
campaign described in Til, but the low quality its H-band 
spectrum prohibited a reliable Mbh estimate, and thus it was 
not included in Tl 1. 

So far, 33 sources of the total of 41 z — 4.8 sources have 
been observed with the SPIRE instrument (iGriffin et al.l2010l) 
on-board Herschel, providing images at 250, 350, and 500 
f(m, corresponding to the rest-frame far infrared (FIR) wave- 
length range of 43-86 jim. All SPIRE observations were 
made in the small-map mode, most suitable for point sources. 
The data reduction process starts with the level 0.5 product of 
the SPIRE pipeline. We then apply the standard tools, which 
are provided by the Herschel Science Centre (HSC) via the 
HIPE software (Ott 2010, version 7.3), and using the version 
7.0 of calibration files. Since all of our sources appear as point 
sources in the images, we follow the guidelines of the HSC 
and apply a peak fitting method in order to measure the to- 
tal flux. This is done by fitting a two dimensional fixed-width 
Gaussian function to the image and taking the value at its peak 
to be the flux of the source. 

In seven sources, the peak emission is shifted by 10-30" 
from the optical location of the AGN. The low resolution 
SPIRE images do not allow proper source separation in such 
cases and deep observations at shorter wavelengths are needed 



to confirm these detections. A Spitzer campaign is underway 
to secure such observations, therefore we defer the analysis 
of these sources to a future publication. Of the remaining 26 
sources, five were detected at above a 3-<7 significance level 
in at least two SPIRE bands, and twenty were not detected at 
all SPIRE bands. One additional source (J 1306+0236) was 
detected above a 3-CT level only in the 350 fim band. In this 
Letter we focus on two groups of sources either detected at 
least in two bands or undetected in all three bands. Therefore, 
we defer the analysis of this source to a future publication and 
hereafter refer only to the remaining 25 sources. 

There are three sources of uncertainty related to the SPIRE 
images (see the SPIRE observers' manuafl). The first uncer- 
tainty is in the fitted value and includes instrument and con- 
fusion noise. The integration time of each observation (either 
222 or 296 seconds) is long enough to minimize the instru- 
ment noise and obtain an image that is dominated by the ex- 
tragalactic confusion noise. The confusion noise is estimated 
to be 5.8, 6 3 and 6.8 mJy/beam at 250, 350 and 500 |iin, re- 
spectively dNguven et alj|2010h . The pixelization correction 
of the images introduces an uncertainty of about 2-3% of the 
flux density. An additional uncertainty is associated with the 
calibration process and is about 7% of the flux density. These 
uncertainties are added in quadrature and listed together with 
the measured fluxes in table Q] 

2.2. SED Fitting and Stacking Analysis 

The 350 lira images of the five detected sources are shown 
in Figure Q] and their spectral energy distributions (SEDs) are 
shown in Figure [2] We also show the SDSS, H-, and K-band 
(where available) spectroscopy from Tl 1 and the m ean intrin- 
sic mid infrared (MIR) AGN SED from lMor & Netzerl(l2012l) . 
The normalization of the MIR SED uses the known Lagn and 
assumes a total covering factor by hot and warm dust of 0.5. 
As discussed in Mor & Netzerj (120121) . this value is close to 
the upper limit for such luminous AGNs. The AGN contri- 
bution to the emission at wavelengths longer than ^30 fim is 
small and its effect on the measurement of Lsf is negligible. 

To measure Lsf we fit the Herschel data points of each 
of the detected sources using a grid of tem plates that span 
a wid e range of FIR luminosities and SFRs dCharv & Elbazl 
2001). The FIR luminosity is the only free parameter for this 
set of templates. Lsf is calculated by integrating the best fit 
model between 8 and 1000 jim. Comparison with other types 
of models which are often used to estimate FIR luminosities 

4 http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html 
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FIG. 1. — 350 /im images of the detected and stacked sources (see text). 
Each image is 2.5' x 2.5' and the blue cross represents the optical (SDSS) 
location of each of the individually detected sources. 

(e.g., gray body) is beyond the scope of this Letter. The fitting 
procedure uses a standard % 2 minimization to determine the 
best fit template. The uncertainty is calculated using the stan- 
dard confidence levels for % 2 with two degrees of freedom. 

As noted, twenty sources were not detected above 3-(7 level 
at any of the bands. To get an average luminosity of these 
sources, we applied a stacking analysis to all twenty images 
in each band. We first cut each image to a small stamp sym- 
metrically around the center of the pixel in which the optical 
location of the source lies. All stamps have an equal num- 
ber of pixels and are approximately 1' x 1' in size. A stacked 
image is constructed by assigning the images with weights 
according to their respective exposure times, and averaging 
the images pixel by pixel. Since the dominant source of un- 
certainty is the extragalactic confusion noise, the slightly dif- 
ferent exposure times have a negligible effect. The stacking 
procedure revealed a statistically significant signal at all three 
bands (hereafter "stacked source") that represents the average 
flux of the individually undetected sources. We then measure 
the flux of the stacked source by fitting a two dimensional 
Gaussian function as described above. 

The averaging of the images assumes that all objects are 
located at the center. Any contribution from neighboring 
sources would be significantly reduced due to the fact that 
such contributions are expected to be randomly distributed in 
the images. The average values can be biased if the (unde- 
tected) sources have very different fluxes. In such a case, few 
sources that are just below the confusion noise limit might 
skew the result towards a higher average flux. To overcome 
this we used a bootstrap approach to estimate the true value of 
the flux of the stacked source and its uncertainty. Out of the 
20 images we choose 10000 random multisets of length 20. In 
each multiset, images may appear more than once. We then 
stacked each multiset of images and measured the flux of the 
emergent stacked source. The probability distribution func- 
tions (PDFs) of the bootstrap procedure, in all three SPIRE 
bands, are shown in Figure [3] The final value of the flux in 
each band is taken to be the maximum likelihood value of 
the corresponding distribution. Its uncertainty is estimated by 
measuring the 16th and 84th percentiles, which are assumed 
to represent the 1-c error. The measured fluxes of the stacked 
source are listed in Table Q] To measure the Lsf of the stacked 
source we follow the same fitting procedure described above. 
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FIG. 2. — Spectral energy distributions (SEDs) of the 5 individually de- 
tected sources. The Herschel/SPIKE detections are shown as circles in 
each panel, and the best fit luminous SF templates are shown as solid blue 
lines. The black lines at the short wavelength end represent SDSS spec- 
tra supplemented with H- and K-band spectroscopy (where available) from 
Trakhtenbrot et al. (2011). Red dashed lines are mean intrinsic MIR AGN 
SEDs from Mor & Netzer (2012) (see text). 



3. RESULTS AND DISCUSSION 

Table [T]lists Lagn (taken from Til), Lsf. an d SFRs for the 
detected and stacked sources. SFRs were calculated from Lsf 
assuming a Salpeter (1955) initial mass function (IMF) and 
using SFR= 1.73 x 10~ 10 Ls F /L w (|Kejmicjit3[1998j). A top- 
heavy Kroupa IMF (Kroupa 2001) would result in SFRs that 
are lower by a factor of ~ 1 .6. 

We have collected from the literature FIR and sub- 
mm measurements of other luminous AGNs that imply 
Lsf > 5 x 10 12 L n . Eight source s at z ^4 have a single 850 
fim detection (Isaaketal. 2002). Four sources have been 
reported by iPriddev et al.l (120031) . Two of these at z ~5 
and z ~6 have a single 850 fim detection and 2 others, at 
similar redshifts, have both 850 /im and 250 GHz detections. 
Eight additional s ources at z ~ 5 .8-6.2 were detected at 250 
GHz dWang et alJl2008l l20lll) . Finally, iLeipski et all d20l0h 
reported multi-band measurements of two sources at z ~4.7 
and 6.4. We used the fluxes reported in these papers and 
applied our fitting method to obtain Lsf- Fitting SF templates 
to one or two data points is problematic in two ways. First, 
a single data point fit does not allow the calculation of a 
meaningful confidence limit. Second, the 850 fim and 250 
GHz bands translate to about 170 jJ.m rest-frame wavelength 
at the reported redshifts. This wavelength is far from the 
peak wavelength emission of cool dust. To estimate the range 
of possible Lsf for the sources with a single data point, we 
fitted the single-band measurement with a gray body with 
emissivity index of j3 = 1.5 and two temperatures, 40 and 
60 K. These temperatures represent the uncertainty in Lsf 
for these sources. Laqn for these sources were calculated 
using the continuum flux density at 1450 A and the relation 
log^AGN/lO^ergs- 1 ) = 0.941og(L 145 o/10 46 ergs- 1 ) + 
0.53, where L1450 = AL^(1450A). This relation is based on 
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FIG. 3. — Probability distribution functions of the fluxes of the stacked 
sources (thin solid blue lines). The vertical solid lines mark the peaks (the 
most likely value) of the distributions. Inner and outer dashed lines represent 
the different percentiles corresponding to l-o" and 3-o~ uncertainties, respec- 
tively. 



several samples of high-z high L A qn sources dShemmer et al.l 
l20^lNelzerirani2007aI andTll). 

Fig. |4] summarizes the main results of this study on the 
Lsf vs. Lagn diagram. It shows a collection of several 
samples of low and interm ediate luminosity AGN-dominated 
sources from Netzer (2009), the high luminosity high redshift 
sources collected from the literature, and our z — 4.8 sources. 
The AGN-dominated correlation line and a 1 : 1 line are also 
shown. 

The two main results of our work are best explained by con- 
sidering the data shown in the insert of Fig.|4] First, the SFRs 
of the five individually detected sources are comparable to the 
most FIR luminous objects known, including SMGs. Thus, 
the most intense SF phase in AGN hosts is not necessarily 
associated with an obscured SMBH phase. Second, the de- 
tected and stacked sources occupy two distinct locations on 
the Lagn-^sf plane. The mean SFR of the detected sources 
is higher than the SFR of the stacked source by ab out a factor 
of 5. T his result is consistent with the findings of Wang et al. 
d2011l) for a sample of z^6 quasars. All the sources in our 
sample have similar optical-UV spectra, span a narrow range 
of Mbh values, similar high values of Lagn, and almost iden- 
tical redshifts (see Table Qjl. Thus, the AGN properties do not 
explain the different SF properties. 

We consider three scenarios to explain the observed differ- 
ent SFRs in the two groups. 

1. The different locations represent different evolution- 
ary routes. The detected sources are in major mergers 
with extremely high SFRs and the undetected ones go 
through a calmer secular evolution with lower SFRs. 
The similar AGN properties and different SFRs of the 
sources in both groups imply different relative growth 
rates of the SMBH and stellar mass in each group. As- 
suming that the ratios of growth rates remain the same 
for some time (before or after z=4.8), this will lead to 



FIG. 4. — AGN luminosity, Lagn, vs. star formation luminosity, Lsf for 
several samples. The five z ~ 4.8 sources detected using the Herschel/SPIKE 
are shown as blue squares. The value of the stacked source is represented 
by the blue empty square. The horizontal error bar represents the range 
of Lagn for the group of undetected sources, and the vertical error bar is 
set by the 3-0" confidence level of the PDFs in Fig [5] Dots represent data 
from several samples of type-I AGN at low and high redshifts (black dots) 
and type-II AGN (smaller red dots), see text. The main correlation for 
AGN dominated sources is shown as a dashed line and can be expressed as 
Lsf — 10 43 (LAGN/(10 43 ergs~')) 07 . A 1:1 ratio is shown as a solid line for 
comparison. Open gray circles represent observations collected from the 
literature. Their Lsf determination is based, in most cases, on a single sub- 
mm point (see text). The insert shows Lagn vs. SFR for the new z — 4.8 data. 

different M*/Mbh at the end of their evolution. Since 
all the sources are probably on their way to become 
the most massive black holes in the universe (Til), this 
scenario is not supported by observations in the local 
universe. 

2. All the sources are in merging systems where an initial 
burst of SF is followed by accretion onto one or two 
SMBHs. The different SFRs may be explained by dif- 
ferent conditions in the mergers (e.g. galaxy size and 
gas content). If this is correct, the SFRs of the sources 
should have been more evenly distributed across the 
SFR range. Furthermore, in this case the gas supply 
to the SMBH is only loosely connected with the SF in 
the host. 

3. All the sources reside in merging systems of similar 
type and the different SFRs are due to different stages 
of the merger process. The process can continue in two 
different ways. The first is that both SFR and L/Leaa 
keep their observed value for a certain time, i.e. Lsf 
remains roughly constant while Lagn increases with 
time. In terms of the Lagn-^sf diagram, it means 
that sources start their fastest SF growth phase at some 
point in the diag ram and travel hor izontally, similar to 
the suggestion of Shao et al. (2010). However, if our 5 
individually detected sources were to follow this path, 
they would reach the AGN-dominated correlation line, 
at lower redshift, with La_gn> 10 49 ergs~'. Such lumi- 



nous AGN are extremely rare at z ~2, and have never 
been observed at z >3 but in our sample they represent 
25% of the population. The second possibility is that 
the five detected sources are at the peak of the SF phase 
and the other 20 towards its end. Given the similar red- 
shift of the sources, the SF must have been reduced sig- 
nificantly over a short period of time, while their Lagn 
remains at its peak. The redshift range of the sam- 
ple (4.65<z<4.92) implies a timescale of —100 Myr 
for this process. In terms of the Laqn-Lsf plane, this 
means that the 20 sources with lower Lsf had higher 
SFR in the past but traveled vertically down towards 
the AGN-dominated correlation line. A possible way 
to rapidly quench SF is by AGN feedback. In this sce- 
nario, which we consider more probable, feedback re- 
duced SF in the hosts of the 20 lower- Lsf objects but is 
not yet operating efficiently in the other five. If this sce- 
nario is true the AGN activity should be significantly 
reduced shortly after the quenching of the SF. This is 
consistent with the estimate of the AGN duty cycle by 
Til. 



To conclude, our optically-selected, flux-limited sample, 
which consists of sources that span a narrow range of AGN 
properties, provides an excellent test-bed of various SMBH 
and galaxy evolution scenarios atz ~ 4.8. The extremely high 
SFRs found in five sources provide strong evidence for a 
merger process in these systems. The clear separation into 
two groups with SFRs that differ by a large factor provides, 
perhaps, indications for AGN feedback. 
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